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We study numerically stabilization against ionization of a fully correlated two-electron model 
atom in an intense laser pulse. We concentrate on two frequency regimes: very high frequency, where 
the photon energy exceeds both, the ionization potential of the outer and the inner electron, and 
an intermediate frequency where, from a "single active electron" -point of view the outer electron is 
expected to stabilize but the inner one is not. Our results reveal that correlation reduces stabilization 
when compared to results from single active electron-calculations. However, despite this destabilizing 
effect of electron correlation we still observe a decreasing ionization probability within a certain 
intensity domain in the high-frequency case. We compare our results from the fully correlated 
simulations with those from simpler, approximate models. This is useful for future work on "real" 
more-than-one electron atoms, not yet accessible to numerical ah initio methods. 



I. INTRODUCTION 

The advent of high intensity lasers led to an increasing interest in non-perturbative studies of atomic systems 
interacting with intense laser light (see, e.g., Q for a review). One of the most frequently revisited topics during the 
last fifteen years was stabilization of atoms (or ions) against ionization in intense laser light, i.e., for increasing laser 
intensity the ionization rate decreases. This kind of stabilization was predicted by Gersten and Mittleman already in 
1975 Experimentally, stabilization of highly excited atoms has been reported ||^ whereas measuring stabilization 
of atoms initially in the ground state is hard to achieve. This is due to the fact that, in order to see stabilization, 
the laser photon energy has to exceed the ionization potential. Unfortunately, there are not yet high intensity lasers 
available delivering such energetic photons. Therefore most of the studies in this field are of analytical or numerical 
nature: "high-frequency theory" Q, Floquet calculations [^ E] , the numerical treatment of ID model atoms, quantum 
[0-0 and classica l [pd| , as well as in two-color laser fields 2D atoms in arbitrary polarized laser light jl^, and 
full 3D hydrogen [|14|. Of particular interest is whether the atom survives in a "real" laser pulse up to intensities 
where stabilization sets in, or whether it already ionizes almost 100% during the rise time of the pulse ||l^. In other 
words: is the atom able to pass through the "death valley" of ionization before arriving at the "magic mountain" of 
stabilization? There are also several papers where the authors came to the conclusion that stabilization does not exist 
at all (see [^0, and references therein). 

In this paper we focus on how the electron correlation in a two-electron model atom affects the probability for 
stabilization, i.e., the probability that the model atom remains neutral after the pulse has passed by. For two 
frequency regimes we compare the results from the fully correlated calculation with approximate models like "single 
active electron" or time-dependent density functional theory. The purpose of these studies is, on one hand, to gain a 
qualitative picture of the stabilization mechanism in a more-than-one electron atom, and, on the other hand, testing 
approximate methods before applying them to 3D many-electron atoms where accurate, full ah initio studies are not 
possible with current days computers. To our knowledge only a few other numerical studies of correlated two-electron 
systems in the stabilization regime are reported in the literature so far |l8|-po| . 

II. THE MODEL ATOM 

We study a model helium atom where both electrons are allowed to move in one dimension only, but with the 
electron-electron correlation fully taken into account. This leads to a two-dimensional time-dependent Schrodinger 
equation (TDSE) 

4l*W) = H(t)|*(t)) (1) 
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with the Hamiltonian 



m = i(p, + Ait)r + J(P2 + Ait)r + \ (2) 



Here, the laser pulse is coupled to the atom in dipole approximation through the vector potential A[t). Xi and Pj 
(i = 1,2) are the electrons' coordinates and canonical momenta, respectively. We use atomic units (a.u.) throughout 
this paper. The regularization parameter e was chosen 0.49 which yielded, on our numerical grid, ionization potentials 
similar to real helium (0.9 a.u. for the first electron, and 2 a.u. for the second one). 

The electric field E{t) = —dtA{t) was a trapezoidal pulse with a rising edge over 5 optical cycles, 5 cycles of constant 
amplitude E, and a down-ramp over, again, 5 cycles. We started all our simulations with the field-free ground state 
|4'(0)). The wavefunction "^{xi, X2, t) = {xiX2\'^{t)) was propagated in time using an unconditionally stable, explicit 
"grid hopping" algorithm [El|. Non-vanishing probability amplitude 'i>{xi,X2) near the grid boundary was removed 
through an imaginary potential. The numerical grid was always several times (at least 10 times) larger than the 
excursion length 

a^\E/c.'\ = \A/u;\ (3) 

of a classical electron oscillating in the laser field of frequency ui and electric field amplitude E (vector potential 
amplitude A). During time propagation we monitored the amount of probability density j^E'p inside a box xi,X2 G 
[— 5,-f5]. After the pulse is over, the density inside this box can be interpreted as the "survival" probability of the 
helium atom to remain neutral [ p2| . 

To analyze the results obtained with this fully correlated model atom we compare with several simplified models. 
Among those, the "single active electron" (SAE) approximation is the simplest one. There, one assumes that an 
inner and an outer electron respond independently to the laser field. The inner electron "feels" the bare nucleus 
{Zi — 2, hydrogen- like) . The outer one sees an effective nuclear charge, to be adjusted in such a way that the correct 
ionization potential (0.9 a.u.) is obtained. In our numerical model this was the case for Z^ = 1.1. Thus, in the SAE 
approximation, we solved two independent TDSEs with no dynamic correlation at all, 

idt^.ix, t) = (\{-id, + A{t)f - -^=\ *,(x, t), (4) 
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idt^oix. t) - ( i(-i5, + A{t)f - -j^= \ ^>o{x, t). (5) 
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In order to incorporate correlation in a first step one can introduce a Hartree-type potential into the Hamiltonian 
for the inner electron, 

ia...(.,t) ^ (i(-ia. + ^(.))^ - + / J^=Ld.') ..(.,.), (6) 

idt^o{x,t) = (h-id. + A{t)f - -^^=\ *o(x,t). (7) 
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In this approximation, the inner electron feels the bare nuclear potential and the outer electron. Therefore, we call 
this model "inner sees outer" (ISO) approximation. It was utilized in Ref. to study non-sequential ionization 
(NSI). In the ground state, the Hartree-potential leads to a screening of the bare nuclear charge. Thus, energetically 
the two electrons are almost equivalent in the beginning, though we labelled them "inner" and "outer" in Eqs. (|^), 
(^. However, during the interaction with the laser field one of the electrons might become the outer one. We will 
also consider the opposite point of view where the outer electron sees the inner one ("outer sees inner", OSI). In this 
case we have to deal with the system of TDSEs 

\dt^,{x,t) = {\{~\d,.+ A(t)f - ^^=\-^,{x,t), (8) 



,2^ ^' ' " ^. 



with Z = 1. 
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Finally, another way to study our model system is to apply time-dependent density functional theory (TDDFT) 
(see ||2^ ] for an overview) in local density approximation, leading to the (nonlinear) TDSE for one Kohn-Sham orbital 
$(x,t), 

idMx,t)=(h-id, + A{t)r-^^+ I -^M^^=Ax']^{x,t). (10) 
\^ Wx^+e J ^J{x ~ x'Y + e J 

The total electron probability density is given by n{x,t) = 2|$(x,t)p. Since, strictly speaking, ^{x,t) cannot be 
interpreted as a physically meaningful single electron orbital it is not easy to deduce single electron quantities (such 
as single ionization for instance). On the other hand, the somewhat arbitrary distinction between an inner and an 
outer electron is avoided in TDDFT. 



III. HIGH-FREQUENCY RESULTS 
A. Single active electron-approximation 



In this Section we want to compare the results from the fully correlated model atom with those from the correspond- 
ing SAE calculations. First, we focus on the high-frequency regime where the laser frequency exceeds both ionization 
potentials, lo = tt > £i > So- From an SAE-point of view we expect both electrons to stabilize, especially the outer 
one since the frequency is 3.5 times larger than the ionization potential So = 0.9. In Fig. |l| the amount of probability 
density inside the box xi,X2 G [—5, -1-5] (PDIB) vs. time for a = 0.5, 1.0, 1.5 and 2.0 is shown for the inner and the 
outer electron. First we observe the expected result that the outer electron is more stabilized than the inner one: the 
amount of PDIB after the pulse has gone is greater for the outer electron. However, qualitatively the set of curves 
are quite similar. For the two higher d-values (1.5, drawn dashed, and 2.0, drawn dashed-dotted) one observes that 
the curves bend sharply after the up-ramping at t = 5 cycles, i.e., for both electrons ionization is much slower during 
the constant, intense part of the pulse. In fact, ionization happens almost exclusively during the rampings. We also 
note that the slight decrease of the PDIB during the constant part of the a = 1.5 and a = 2.0-pulses is linear in time, 
in contrast to tunneling ionization where we have an exponential dependence (in the case of a constant ionization 
rate). Approximately two cycles after the down-ramping at t = 10 cycles, ionization starts to increase again. Finally, 
after the pulse is over (at t = 15 cycles) the amount of PDIB remains stationary. Secondly, we observe that there is 
obviously no monotonous behavior of stabilization with increasing intensity (or a). Ionization is higher for a = 1.0 
(dotted curves) than for a = 0.5 (solid curves). For a = 1.5 (dashed curves) ionization starts to decrease, i.e., we are 
entering the stabilization domain at that point. Therefore the so-called "death valley" seems to be located around 
a = 1 for both electrons in our model atom. 

In Fig. H the stabilization probability for both electrons is shown vs. the excursion a. The quiver amplitude a is 
related to the laser intensity I — through a = The stabilization probability of the inner electron exhibits 

an oscillatory behavior. The "death valley" is located at d « 1, followed by a maximum at d w 4. For higher intensity 
ionization increases again up to a stabilization minimum around a ~ 7. Stabilization of the inner electron recovers 
till the next maximum at d w 9. The second maximum is below the first. Thus we observe an overall decrease of 
stabilization with increasing intensity. This is even more pronounced in the stabilization probability for the outer 
electron where the oscillations are less visible. The "death valley" for both electrons is located at d « 1 while the 
maxima are at different positions. The oscillatory character of the stabilization probability in ID systems has been 
observed by other authors as well p5|-|27t . In contrast to our results an overall increase of stabilization with increasing 
intensity was found in |^,^. This might be due to the fact that we are looking at the ionization probability after 
the pulse is over while in analytical papers often the ionization rate is discussed. In the former, ionization during the 
up and down-ramps is taken into account while in the latter it is commonly not. 

The probability for our He model atom to survive as neutral He after the pulse is over is, in SAE approximation, 
simply the product of the probabilities for each electron to remain bound. In Fig. the corresponding curve is 
indicated by i-o. The result from the fully correlated system is also shown (drawn dotted, indicated with 'corr'). We 
infer that, especially for d < 5 the stabilization probability is strongly overestimated by the SAE treatment. We could 
argue that if the system stabilizes it will probably stabilize in such a way that the correlation energy is minimized. In 
that case it sounds more reasonable to take the square of the SAE stabilization probability for the inner electron. In 
what follows we will refer to this viewpoint as "independent electron" (IE) model since it follows from crossing out the 
correlation term in the Hamiltonian (|). The result is included in Fig. ^, labelled i-i. The IE-curve seems to oscillate 
around the fully correlated one, especially for d > 6. From this result we conclude that, compared to the IE-model 
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with two equivalent inner electrons, electron correlation washes out oscillations in the stabilization probability and, 
therefore, can stabilize as well as destabilize, depending on the intensity (for a given pulse shape and frequency). 

To discuss that further we look at the time-averaged Kramers-Henneberger potential (TAKHP), i.e., we transform 
to the frame of reference where the quivering electron is at rest but the nuclear potential oscillates, and average over 
one cycle, 

yKH{xi,X2) = Y ^ I ' At+ , ^ =. (11) 

^2^io 0-49 ^{xi - X2Y + 0.49 

For sufficiently high frequencies this is the leading term in a perturbation series in Q. In the correlation term 
no a{t) appears since the interparticle distance is not affected by the KH transformation. We calculated numerically 
the TAKHP. The result is shown in Fig. ^ for a{t) = asmwt with d = 5. For comparison the TAKHP with the 
correlation term neglected is also shown (corresponding to the IE model). With correlation, there are two minima 
near xi = a^X2 — ~ol and x\ = —a,X2 = a whereas without correlation there arc two more, energetically equivalent 
minima at xi = X2 = a and xi — X2 — ~a. However, if we assume that the fully correlated system manages it 
somehow to occupy the ground state of V^°^^, the correlation energy will be small (for not too small d) since the 
interparticle distance is 2d. The higher d the lower the correlation energy. We believe that this is the physical reason 
that, for increasing d, the agreement of the IE results with the fully correlated ones becomes quite good (although the 
latter do not exhibit an oscillating stabilization probability). Our viewpoint is further supported by examining the 
probability density of the fully correlated system during the pulse. In Fig. || [^'(xi, a;2)P is shown for w = tt, d = 4.0 
at t = 7.5 cycles, i.e., in the middle of the constant part of the trapezoidal pulse. We clearly observe dichotomy, 
i.e., two probability density peaks at the classical turning points, well known from one-electron systems Due to 
electron correlation we do not observe four peaks. Instead the peaks at xi ^ X2 = ±4 are suppressed, in accordance 
with our discussion of the TAKHPs in Fig. S. Therefore the correlation energy is rather small since the distance 

is V8d 



between the two peaks in the a;ia;2-plane is v8d. In the work by Mittleman |2§f| such multi-electron "dichotomized" 
bound states are calculated. 



B. Time-dependent density functional theory 



In Fig. Hour results from the TDDFT calculations are presented. Although the Kohn-Sham orbital ^{x,t) is an 
auxiliary entity that has, in a rigorous sense, no physical meaning, we take it as an approximation to a single electron 
orbital. If we do this, \^{x, t)p • \^{x, t)p, integrated over the region —5 < a; < 5 after the pulse is over, is our TDDFT 
stabilization probability. We see that for d < 1.5 the agreement between TDDFT and correct result is very good. 
The difference between TDDFT and IE (indicated by i-i again) is a direct measure of correlation effects since both 
models differ by the Hartree-term in the Hamiltonian only. Up to d « 5.5 electron correlation suppresses stabilization 
compared to the IE approximation. In that region TDDFT agrees better with the full result. As mentioned before, for 
higher d the IE curve oscillates around the correct result and therefore it comes occasionly to a very good agreement 
with the exact result. Also the TDDFT result agrees very well with the fully correlated curve for d > 7. In summary 
we can say that the TDDFT result is in good agreement with the exact, fully correlated stabilization probability. Both 
have their maximum around d ~ 4 and the "death valley" is also at the right position. For higher d the agreement 
seems to become even better. 



C. "Inner sees outer" and "outer sees inner" -approximation 

In order to explain non-sequential ionization (NSI; see, e.g., Ref. [|l| for an overview) it is essential to incorporate 
electron correlation (see |^ for a very recent paper, and references therein). For that purpose Watson et al. added 
a Hartree-type potential to the TDSE for the inner electron (see Eq. ||, ISO). By doing this the double ionization yield 
is greatly enhanced, in accordance with experimental results [ |30| ] . The question we address in this Section is whether 
this method is applicable to stabilization as well. We will also study the opposite procedure, i.e., where the outer 
electron feels the inner one (see Eqs. ^ and ^ OSI). From the discussions on the SAE approximation above we can 
expect that ISO will probably not agree very well with the exact results since the assumption that the outer electron 
sees a static, effective nuclear charge is not valid. In Fig. |6| we see that, for low d ISO (A) behaves like SAE (the 
i-o-curve) while OSI (O) is similar to IE (the curve indicated by i-i). In ISO approximation for d > 3 the electron 
correlation obviously causes strong ionization, compared to the SAE result. Especially during the down-ramping, 
when probability density of the outer electron moves from the turning points ±d back toward the nucleus, ionization 
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of the inner electron is enhanced. For a > 4.5 the ISO curve even drops below the exact result (indicated by 'corr'). 
In OSI approximation the stabilization probability is also underestimated for a > 4.5. In summary we can say, 
that for a < 2, OSI is in very good agreement with the correct result while ISO is not, due to the inappropriate 
assumption of the outer electron feeling just a static effective nuclear charge. Howeve r, for higher a ISO and OSI 
tend to underestimate the stabilization probability while TDDFT does not (see Section [II B and Fig. 0). 



IV. INTERMEDIATE FREQUENCY RESULTS 

In this Section we discuss the stabilization probability in the intermediate frequency regime So < oj < £i where, 
according a single active electron point-of-view, the outer electron should stabilize while for the inner one ionization is 
more likely. In Fig. |^ we compare the result from the fully correlated calculation with those from the SAE treatment. 
In SAE approximation, the outer electron is more stable than the inner one in the region 1.5 < a < 8.5. For the 
inner electron no clear stabilization maximum is visible. For the outer electron the maximum is at d ~ 6, i.e., it is 
shifted toward higher a compared to the high-frequency case. Both, i-o and i-i underestimates ionization, especially 
for low a in the "death valley" -region. Electron correlation obviously enhances ionization. For lower frequencies this 
is the well-known effect of NSI ( and references therein). Although in the fully correlated result we observe a 
stabilization probability maximum around a « 7.5 the absolute value is below 0.04, and, in our opinion, it makes no 
sense to talk about real "stabilization" in that case. As in the high-frequency-case we observe an overall decrease of 
the stabilization probability for very high d- values. 

In Fig. H we compare the result from the fully correlated model atom with the corresponding ones from the ISO, 
OSI, and TDDFT runs where electron correlation is included approximately. Let us first focus on the d-region "left 
from the death valley" , i.e., d < 1. There we observe that ISO is nearest to the correct result while OSI and TDDFT 
underestimates ionization. This is quite understandable within the present knowledge of how NSI works: the inner 
electron needs to interact with laser field and the outer electron in order to become free. Obviously this is best 
accounted for in the ISO approximation. For most d TDDFT lies between the ISO and the OSI result. This is also 
quite clear since in TDDFT both correlated electrons are treated on an equal footing (one Kohn-Sham-orbital only) 
whereas in ISO (OSI) the inner (outer) electron feels the outer (inner) partner through Coulomb correlation, but not 
vice verse. However, all these approximations still overestimate the stabilization probability, at least in the interesting 
d-regime where the stabilization probability rises at all (i.e., for 2 < d < 7.5). 

To summarize this Section we can say that in order to achieve stabilization of our two-electron model atom it is 
necessary to choose a laser frequency that exceeds all ionization potentials. For an intermediate frequency the outer 
electron cannot stabilize owing to correlation. The SAE picture is not appropriate and even ISO, OSI, or TDDFT 
where electron correlation is included approximately fail. 



V. DISCUSSION AND SUMMARY 



In this paper we studied how the electron correlation in a two-electron model atom affects the probability for 
stabilization. We found clear stabilization only for frequencies that exceed both ionization potentials. Although 
for the intermediate frequency we did not find a monotonous increase of the ionization probability with increasing 
intensity we prefer not calling this effect stabilization since, on an absolute scale, its probability was very small. In all 
cases electron correlation reduced the stabilization probability compared to the SAE picture. In the high-frequency 
case the two electrons behave more like two independent inner electrons. Similar results were obtained by Grobe and 
Eberly [^sj for a H^ model-ion. Lewenstein et al. [|9| performed classical calculations for a model-atom similar to 
ours which also showed that "dichotomized" two-electron states are dynamically accessible. 

The agreement of the exact numerical result with TDDFT in the high-frequency case was quite good while in 
the intermediate frequency regime stabilization was overestimated by all approximate techniques (ISO, OSI, and 
TDDFT). 

It is well-known that a slow time-scale in the stabilization dynamics is introduced owing to floppings between states 
in the time-averaged Kramers-Henneberger potential It can be easily imagined that these slow floppings are 

affected by electron correlation because, e.g., the merging of the two dichotomous peaks into a single one is suppressed 
then. But even without correlation the results for the stabilization probability are quite sensitive to rise time and 
pulse duration of the laser field since it strongly depends on which Kramers-Henneberger states are mainly occupied 
at the time instant when the laser pulse ends. To avoid these additional complications in the interpretation of the 
numerical results we chose a rather short laser pulse duration so that low-frequency Rabi-floppings do not play a 
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role. Therefore, in the high-frequency studies we just observed the two dichotomous peaks building up (as depicted 
in Fig. ^) but no peak-merging during the constant part of our trapezoidal laser pulse. 

Finally, we would like to comment on the reduced dimensionality of our two-electron model atom. We also performed 
calculations with "real", i.e., three-dimensional (3D) hydrogen- like ions in the stabilization regime. It seems to be the 
case that in 3D stabilization is less pronounced. Moreover, the oscillatory character is less visible, i.e., we observe a 
single stabilization maximum followed by a rather monotonous increase of the ionization probability. The difference 
of ID models and 3D hydrogen was also studied in Ref. |2^. The effect of electron correlation in 3D stabilization will 
be the subject of a future paper |31|. 
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FIG. 1. The amount of probability density inside the box Xi,X2 G [—5, +5] vs. time for a — 0.5, 1.0, 1.5 and 2.0 (drawn solid, 
dotted, dashed, and dashed-dotted, respectively) for the inner and the outer electron in "single active electron" -approximation. 
The laser pulse of frequency oj = tt was ramped-up linearly (in field) over 5 cycles, held constant for another 5 cycles before 
the linear down-ramp between t = 10 and t = 15 eye. As expected, the outer electron is more stabilized. For the two higher 
d- values we clearly see that ionization is weak during the intense, constant part of the pulse. 



FIG. 2. Stabilization probability for inner (i) and outer (o) electron in SAE-approximation for ui = n. The fully correlated 
treatment leads to the stabilization probability (i.e., the probability for the He model atom to survive as neutral He) drawn 
dotted and indicated with 'corr'. From a SAE viewpoint one would expect the product of the curves i and o, indicated with 
'i-o'. Obviously, in SAE approximation the stabilization probability is overestimated. The so-called "death valley" (dv) for 
small a is located around a ~ 1 for both, SAE and fully correlated results. The i-i-curve results from an independent electron 
model with two inner electrons. This curve seems to oscillate around the correct curve 'corr'. 



FIG. 3. The time-averaged Kramers-Henneberger potential with (V^jj*^) and without (Vkh) electron correlation included 
for Q = 5. With correlation, there are two minima near xi = &,X2 — —a and xi = —a,X2 = & whereas without correlation 
there are two more, energetically equivalent minima at a;i = X2 = d and xi = X2 — ~a. However, once the fully correlated 
system is in the ground state corresponding to V^°j!j^ , the correlation energy is small (for not too small d) since the interparticle 
distance is 2d. 



FIG. 4. Probability density |^'(xi, X2)\'^ for = tt, d = 4.0 at t = 7.5 cycles, i.e., in the middle of the constant part of the 
trapezoidal pulse. We clearly observe dichotomy, i.e., two probability density peaks around xi — ^4:,X2 ~ ±4. Therefore the 
correlation energy is rather small since the distance between the two peaks in the a;ia;2-plane is y/Sa. 



FIG. 5. Stabilization probability in the high-frequency case {uj = tt) calculated from time-dependent density functional 
theory (TDDFT). The fully correlated results (labelled 'corr') and the independent electron curve (indicated by i-i) are included 
for comparison. The agreement of TDDFT with the exact result is good: "death valley" and stabilization maximum are at the 
same d position. For higher d the agreement is even better. 



FIG. 6. Comparison of "inner sees outer" (ISO, A) and "outer sees inner" (OSI, O) results with the fully correlated ones 
(corr). Also shown are the "single active electron" (SAE) stabilizati on pr obability (i-o) and the "independent electron" (IE) 
prediction (i-i), already presented in Fig. ^ and discussed in Section [II A For low d ISO behaves like SAE while OSI agrees 
with IE. For d > 4.5 the ISO curve drops below the exact result. In OSI approximation the stabilization probability is also 
underestimated for d > 4.5. 



FIG. 7. Stabilization probability for inner (i) and outer (o) electron in SAE-approximation for the intermediate frequency 
u! = n/2. The fully correlated treatment leads to the stabilization probability drawn dotted and indicated with 'corr'. From a 
SAE viewpoint one would expect the product of the curves i and o, indicated with 'i-o'. The i-i-curve results from an independent 
electron model with two inner electrons. All these models underestimate ionization, i.e., they overestimate stabilization (please 
note that the stabilization probability is on a logarithmic scale now). 



FIG. 8. Comparison of "inner sees outer" (ISO, A), "outer sees inner" (OSI, O), and time-dependent density functional 
theory (TDDFT, +) results with the fully correlated ones (corr). Stabilization is overestimated in the ISO, OSI, or TDDFT 
approximation. 
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Fig. 1: Bauer and Ceccherini, "Electron correlation vs. ..." 




Fig. 2: Bauer and Ceccherini, "Electron correlation vs. ..." 




Fig. 3: Bauer and Ceccherini, "Electron correlation vs. ..." 
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Fig. 4: Bauer and Ceccherini, "Electron correlation vs. ..." 




Fig. 5: Bauer and Ceccherini, "Electron correlation vs. ..." 




Fig. 6: Bauer and Ceccherini, "Electron correlation vs. ..." 
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Fig. 7: Bauer and Ceccherini, "Electron correlation vs. ..." 
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Fig. 8: Bauer and Ceccherini, "Electron correlation vs. ..." 



